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Abstract 

We report on the possibihty of detecting a submilUmetre-sized extra dimension by observing 
gravitational waves (GWs) emitted by pointhke objects orbiting a braneworld black hole. Matter 
in the 'visible' universe can generate a discrete spectrum of high frequency GWs with amplitudes 
moderately weaker than the predictions of general relativity (GR), while GW signals generated 
by matter on a 'shadow' brane hidden in the bulk are potentially strong enough to be detected 
using current technology. We know of no other astrophysical phenomena that produces GWs 
with a similar spectrum, which stresses the need to develop detectors capable of measuring this 
high-frequency signature of large extra dimensions. 
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Motivation String theory-inspired braneworld models P, 0] envisage our universe as a 
4D membrane embedded in some higher-dimensional space. Standard Model particles and 
fields are assumed to be confined to the 'brane', while gravity propagates in the higher- 
dimensional 'bulk'. The principal observational features of such models take the form of 
modifications to 4D gravity. For static situations, we expect deviations from Newton's law 
at distances less than the curvature scale of the bulk, ^. Precision laboratory measurements 
yield that £ < 0.1mm Similarly, as we demonstrate here, in the case of dynamic 
gravitational fields we can expect significant higher- dimensional effects at frequencies in 
excess of ~ c/^. Presently, there is a vigorous worldwide effort to build detectors capable 
of observing dynamic gravitational degrees of freedom (i.e., gravitational waves) and hence 
verify one of the last untested predictions of Einstein's general relativity. An intriguing 
question is how the braneworld paradigm may affect what these detectors see, and whether 
or not we can obtain useful constraints on extra dimensions. To address these issues, we 
need to concretely model how GWs are generated in braneworld scenarios. The purpose of 
this work is to predict the spectrum and amplitude of GWs generated by pointlike bodies 
orbiting a braneworld black hole. 

The black string braneworld As in our previous work with R. Maartens 0], we model 
a braneworld black hole as a 5D black string spacetime between a 'visible' brane at ?/ = and 
a 'shadow' brane at y = d. The metric is given by ds^ = e~'^^^^^^dsl^-^^ + dy'^, where ds^^^^ 
is the 4D Schwarzschild line element. We require d/i > 5 to comply with post-Newtonian 
solar system constraints This looks like the Schwarzschild solution on 'our' visible brane, 
with deviations from GR appearing perturbatively. This solution is stable if the mass of the 
string M is large compared to the brane separation: M/Mq > 1.1 x 10~^(£/0.1 mm) e^*^"^^)/^. 
If the string is too light an instability shows up in the spherical perturbations, but not for 
higher multipoles '6',7'|. 

Kaluza-Klein radiation from orbiting particles A small compact object (which 
we model as a delta function) of mass Mp on either brane will orbit the black string in the 
same way as it would a Schwarzschild black hole in 4 dimensions. As in GR, it will emit 
massless spin-2 gravitational radiation, but unlike GR it will also generate fluctuations in the 
brane's position as well 'Kaluza-Klein' (KK) modes. These are 5D massless GWs that have 
momentum along the extra dimension, and so behave like a coupled system of spin-0, spin-1, 
and spin-2 massive flelds on either brane. Owing to the flnite separation of the branes, the 
spectrum of KK masses m„ is discrete. We concentrate on the spherical component of KK 
radiation generated by particles orbiting the string on either the visible or shadow brane. 

As described in the Appendix, spherical KK modes are governed by a pair of one- 
dimensional coupled wave equations sourced by the particle. To illustrate typical waveforms 
seen by distant observers, we numerically integrate these equations in the case when the 
dimensionless KK mass is /i = GMm/hc = 0.5. We consider two types of source trajectory: 
an eccentric periodic orbit (Fig. Q), and a 'fly-by' orbit (Fig. |2)). In the former instance a 
nearly monochromatic steady-state signal is seen far from the string, while in the latter case 
we see a burst of radiation followed by a slowly decaying tail. 

For a given source, it is possible to estimate the amplitude /i„ of each massive mode 
as measured on earth. We show these amplitudes as a function of their frequency /„ in 
Fig. El for a number of different cases. Note that the KK frequencies are bounded from 
below: /„, > /mm ~ 12 GHz (£/0.1 mm)~^ e"*^'^"^^-'/^. The variation of amplitude with fre- 
quency is qualitatively different for /„ greater or less than a critical value: /cnt = c/n'^i ~ 
304 GHz [i/O.l mm)~^. For 'visible' sources located on our brane, the amplitudes are peaked 
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FIG. 1: Radiation from a periodic orbit. The steady-state KK gravitational wave signal 
induced by a particle undergoing a periodic orbit around the black string with = 0.5. The orbit 
{bottom left) has eccentricity 0.5 and angular momentum 3.62 in the notation of Refs. 0, 0|. The 
waveform of radiation falling into the black string is quite different than that of radiation escaping 
to infinity: The infalling signal precisely mimics the orbital profile of the source, while the outgoing 
signal is dominated by monochromatic radiation whose frequency is proportional to the KK mass 
/ = mc'^/27rh. 




FIG. 2: Radiation from a fly-by orbit The massive mode (/i = 0.5) GW signal from a fly-by 
orbit (bottom left) with eccentricity 3.0 and angular momentum ~ 3700. We show the generated 
radiation in a spacetime diagram (center), as seen by an observer close to the string {top left), and 
as seen by a distant observer (right). 
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for /„ ~ /crit- For 'shadow' sources located on the other brane, the hn are both independent 
of frequency and maximised for < /crit- In either case, the amphtudes are bounded from 



above hn < /imax, where 



h A 



Mp\ f r \~ f M 



MqJ \kpcj \MqJ VO.lmm 

{5.0 X lO^^^e"'^'^"^^)/^, visible source, 
9.1 X io-2ig~(d-5^)/2^^ shadow source. ^ ' 

Here, r is the distance to the string and A is the characteristic signal amplitude determined 
from simulations. We find that A ^ 1 for fly-by orbits and A ~ 10~^ for periodic orbits. In 
the former case, /imax is comparable to the expected signal strength in GR. 

Detection scenarios Consider a GW detector whose sensitivity is characterised by 
the spectral noise density S{f). If Hn{t) is the detector's (linear) response to the n^^ KK 
mode, then the signal-to-noise ratio for the total massive mode signal built-up over an ob- 
servation time T is SNR = E„ 2S-\fn) Jq Hl{t) dtY'^ 0. Since the frequency separation 
~ /critC"'^''^ between KK modes is small in most cases, this sum can approximated by an 
integral. This leads to the semi-empirical formula: 

SNR ~ Tp-^/^^Z, periodic orbit, 

biNK~ie' /imax < 5/6 .1/3 n , 

[0.83 TJ^ , fiy-by orbit, 

where Tstr = 2GM/c^ ^ 4.9(M/M0) /xs is the characteristic timescale set by the string 
mass, which is typically much less than the observation time. The pre-factor depends on 
the detector, particle orbit, and source location: 



r = r[5,d, 

where 



.1/2 



(3) 



Q{u) periodic orbit fly-by orbit 
visible source 

shadow source (1 + u'^)^^^'^ + -u^)"-^/^ 

Detectors which maximise F stand the best chance of detecting the KK signal. 

The simplest noise model for a GW detector is one in which the characteristic strain 
sensitivity /^strain = \/^{f) is constant over a band [/o — A/, /o -f- A/], and is otherwise 
infinite. In Fig. El we show the /istrain required of such a detector in order to observe 
periodic-orbit KK radiation as a function of /q. We hold the logarithmic bandwidth of the 
detector constant, which yields that the minimum sensitivity is actually independent of /o 
for /o > /crit- That is, the optimal means of detecting this type of KK radiation is via a 
high-frequency GW detector with /o > f cHt. Several designs for devices approaching this 
zone have been proposed or implemented jl^i Fig. |3] shows how such a detector can be 



used to constrain the fundamental parameters of our model. 

Eqs. (H)) and ^ imply that the detector-independent ratio SNR/F decreases exponentially 
with d/£ for visible matter, but is independent of brane separation for shadow matter. 
Furthermore, as seen in Fig. the KK amplitudes hn generated by shadow matter are not 
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FIG. 3: Massive mode amplitudes for visible and shadow sources. We have taken I = 
0.1mm, the string mass to be 10 M©, the particle mass to be 1.4 Mq, source distance to be Ikpc, 
and A = 1. For comparison, the dimensionless one-year sensitivity curve for Advanced LIGO 
(ALIGO) |10] is shown 11 1. The shaded region shows the minimum /igtrain Hz^^^ required to detect 
the KK signal with SNR > 5 as a function of detector frequency /o, assuming a one-year integration 
time and a periodic source. 



suppressed for /„ < /crit, unlike the visible matter case. These facts imply that a direct 
detection of shadow matter is within the capability of 'low-frequency' devices such as LIGO 
and ALIGO, which are otherwise insensitive to KK radiation from realistic sources on our 
brane. In Fig. we show the types of events detectable by these two interferometers when 
they achieve their respective design sensitivities. 

Discussion Gravitational waves may well prove to be a critical tool in placing limits 
on, or detecting positive signatures of, large extra dimensions. By using a simple model 
of a braneworld black hole, we have shown how brane-confined matter can generate GW 
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FIG. 4: Observations in the extra-dimensional parameter space. Signal-to-noise prediction 
for massive mode GWs generated by a fly-by orbit on the visible brane for a high-frequency detector 
with /o = 10^^ Hz, A/ = 10^^ Hz, and /istrain = 10"^^ Hz"^/^. The particle mass is 1 M© while the 
string corresponds to the 'black hole' at the galactic centre (M ~ 4 x 10^ Mq and r ~ 8kpc); the 
scaling of SNR with Mp and r is given by ^ . Also shown are the number of KK modes in this 
detector's range as a function of ^ and d/i. Assuming that the other parameters can be obtained 
from other observations (e.g. the zero- mode GW signal), simultaneous measurements of SNR and 
N can determine the bulk curvature scale and brane separation uniquely. 



FIG. 5: Prospects for detecting KK radiation from shadow matter particles by LIGO 

and ALIGO. For any black string-shadow particle interaction, we can assign a detectability 
statistic X = A{Mp / MQ){r /'kpc)~^ . The contours indicate events that are detectable with SNR > 
1 (assuming a one-year integration time in the case of periodic sources, and ^ = 0.1 mm). We have 
assumed that the detector sensitivity is entirely photon shot-noise/radiation pressure limited for 
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signals with a characteristic frequency fixed by the curvature scale of the extra dimension — 
around 300 GHz or above — and with amplitude comparable to the predictions of GR. The 
most efficient (and cost-effective) means of observing KK radiation generated on our brane 
is with a high-frequency GW detector optimised for all spins of radiation. The detection of 
such a signal would provide a strong evidence in favour of new physics (such as a large extra 
dimension), since there are no astrophysical sources that can generate a similar spectrum of 
GWs. (Certain models of inflation can generate GWs with / ~ 10^^ Hz fxEi, but not with a 
discrete spectrum.) 

It is remarkable that both the intrinsic strength and shape of the shadow matter GW 
spectrum implies that it is easier to detect than KK radiation from visible matter. In 
braneworld models, the only means of directly observing material hidden in the bulk is via 
gravitational interactions, so the search for KK radiation is one of the few viable techniques 
that can constrain the shadow brane's matter content. Intriguingly, our work suggests that 
real observational constraints on shadow matter are within the grasp of current technology 
such as LIGO. 

We note that the discrete frequencies within the KK spectrum are independent of the 
string and particle masses. That is, all such systems will generate GWs with the same 
frequencies, hence there will likely be a signiflcant (non-cosmological) integrated stochastic 
GW background in this model. Indeed, one would expect such a background to be generated 
in any model incorporating 'ultraviolet' modiflcations to GR at length scales i < 0.1mm. 
Hence, stochastic GW backgrounds might prove to be a useful model-independent observa- 
tional window onto submillimetre-scale exotic physics. 

Finally, we reiterate that all of our results have been calculated using the delta-function 
approximation for brane sources. An important open issue involves the effects of more 
sophisticated source modeling, but this is left for future work. 
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Appendix In the Randall-Sundrum gauge [2], perturbations of the black string metric 
are orthogonal to the extra dimension and given by gap — gap + ^n{y)h^api^^)- The 

linearized Einstein fleld equations yield that the Zn are eig enfunctions of -e-^\y\^^{d'^-4:/i^) 

with discrete eigenvalues m\c^ /ti^, which are the effective masses of the KK gravitons h^ap- 
The n = contribution is massless and hence reproduces ordinary GR. To avoid the Gregory- 
Laflamme instability, we need /i„ = GMnin/hc > 0.4301 for all n > 0. 

The individual components of the spherical part of h^^^ can be derived from a master 
variable ip, which satisfles 

(d^^-dl + V^)^ = S^+Iip, (5a) 
(a^ -dl + V^)ip = 5^, (5b) 

where we have deflned the dimensionless coordinates r = tc^/GM and x = rc^/GM + 
2 In (rc^/2G'M — 1). The tortoise coordinate x maps the event horizon at r = 2GM/c^ 
onto X = — oo. In these wave equations, if is another master variable that governs spherical 
perturbations in the position of the brane on which the matter source resides. The equations 
are coupled by the interaction operator X = I{x, dx), and V^^^ are potentials. 

The source terms and depend on the perturbing brane matter. As in GR, we 
analytically model a small brane particle using a stress-energy tensor with delta-function 
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support along its worldline. In numeric simulations, the delta-functions are replaced with a 
narrow Gaussian profile 17 1. There are a few ambiguities in this regularization scheme. 



but we find that our numeric results far from the string are largely insensitive to the par- 
ticular choices made. 

A detailed analysis leads to the following late-time/distant-observer approximation for 
the KK metric perturbations: 

hal - hne'^"' diag (O, +1, -fr^, -ir^ sin^ 6) 

{1, periodic orbits, , , 

(tcVGM)-5/6, fiy-by orbits, ^ ^ 

where \hn\ = VStiA ^^^t^ j (^iP^) Fn{d/l). ^ is a dimensionless quantity that depends 
on the particle orbit but not on any other parameters; its value must be determined from 
simulations. Fn{d/£) is a complicated expression with the following limiting behaviour: 
When the perturbing matter is on our brane 



Fn{d/i) ^ { '„rf/2^. . 1/2 ' ^ ^ (7a) 



e 



On the other hand, for shadow particles 



Finally, to a good approximation, the KK frequencies are given by ujn = 27r/„ 
f(^ + i) 
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